arXiv:hep-ph/9701309vl 16 Jan 1997 


CERN-TH/96-362 
NORDITA 96/79-P 
January 1997 
|hep-ph/9701309 


Phenomenology of Power Corrections in 
Fragmentation Processes in e+e“ Annihilation 


M. Beneke 

Theory Division, CERN, 

CH-1211 Geneva 23, Switzerland 

V.M. Braun and L. MagneaQ 
NORDITA, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 

Abstract 

We analyse power corrections to longitudinal and transverse fragmentation processes in 
e’''e“ annihilation, based on the assumption of ultraviolet dominance of power correc¬ 
tions. Under this assumption, we determine the dependence of power corrections on the 
scaling variable x from the infrared renormalon asymptotics of leading power coefficient 
functions. Our results suggest that the longitudinal and transverse gluon fragmentation 
coefficient functions receive corrections of order 1/{xQY. The power expansion breaks 
down ai X < A/Q and has to be resummed. This resummation leads to 1/Q corrections 
to the longitudinal and transverse cross section, which cancel for the total cross section. 
We provide a simple parametrization of the x dependence of 1/Q^ corrections to frag¬ 
mentation processes and investigate perturbative corrections to the longitudinal cross 
section in higher orders, in view of a determination of the strong coupling. 


On leave of absence from Universita di Torino, Torino, Italy. 






1 Introduction 


Inclusive single-particle production in collisions, —>■ 7 *, 2'° — H{p) + X, 

tests scaling violations in the time-like region and can be used to measure the strong 
coupling Us iii- Additional insight can be gained from a measurement of the angular 
dependence of the detected hadron P, ||, since, for example, gluon fragmentation enters 
the longitudinal cross section as a leading contribution. The differential cross section 
can be expressed as 


dxd cos 9 


(e+e 


HX) 
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(1 -h cos^ 6 ) ^ sin^ 6 

dx 4 dx 


+ |cos0^(a;,g"). (1) 

We dehned x = 2p ■ q/q'^, where p is the momentum of H and q the intermediate gauge 
boson momentum. denotes the center-of-mass energy squared and 9 the angle 

between the hadron and the beam axis. The angular decomposition in (|^ corresponds 
to the contributions from longitudinal and transverse polarizations of the intermediate 
gauge boson and from 7 * — interference. In the following, we will not be concerned 
with the asymmetric contribution and with quark mass effects. Neglecting quark masses, 
{1/ao) daf^lj^ldx (where cxo is the Born total annihilation cross section) is independent 
of electroweak couplings and the longitudinal cross section is suppressed by as- The 
most precise measurements of fragmentation functions refer to a sum over all charged 
hadrons. In the following, we drop the superscript ‘hf’ when a sum over all hadrons H 
is understood. The conversion from charged hadrons to all hadrons is considered as an 
‘experimental problem’. 

It follows from the factorization properties of perturbative QCD that the ‘structure 
functions’ in (|^ are convolutions of perturbative coefficient functions Cp{x, (P = 

T, L, A) and parton fragmentation functions D f{x,p) {i = q,q,g)^, 

—C"p{z,Q‘^/p^)Df{x/z,p), ( 2 ) 

dx ^ Jx z 


up to corrections suppressed by some power of A/Q, where A is the QCD scale parameter. 
The fragmentation functions have to be determined experimentally, but once measured, 
they can be used to predict fragmentation processes in other hard collisions. Both 
coefficient and fragmentation functions depend on the chosen factorization scheme and 
scale /i. We choose dimensional regularization with MS subtractions. Then the coefficient 
functions are obtained as the partonic cross sections with poles minimally subtracted. 
The fragmentation functions satisfy time-like evolution equations with kernels known to 

ipor i = q (or q) we always imply that all light quark flavours are summed over. In the context of 
power corrections a quark is considered ‘light’ if its mass is smaller than A. Charm and bottom quarks 
must be treated separately, even if the center-of-mass energy is much larger than their masses. 
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next-to-leading order p. The coefficient functions in (|^) have been computed to order 

1 - 

The evolution kernels and perturbative corrections to coefficient functions cause scal¬ 
ing violations that are logarithmic in the center-of-mass energy. As in deep-inelastic 
scattering (DIS), multi-parton correlations, not taken into account in (^, lead to scal¬ 
ing violations that scale as some power of A/Q. Such corrections can nevertheless be 
important, as they can be enhanced in specihc kinematic regions. Such a situation is 
well known for DIS, where for XBj —^ 1 all partons must be highly correlated in or¬ 
der that all momentum can be transferred to a single parton. In this region, the twist 
expansion breaks down (see for example Ref. [0). Unlike DIS, power corrections (in 
analogy with DIS we use ‘higher-twist correction’ synonymously) in fragmentation pro¬ 
cesses have rarely been studied theoretically [Q. The main difference comes from the 
applicability of the operator product expansion to DIS, which allows one to express the 
moments of multi-parton correlation functions in terms of matrix elements of local 
operators. Such a classification seems to be more difficult for fragmentation, so that 
even the power behaviour of higher-twist corrections is not well established. Although 
the collinear (light-cone) expansion for fragmentation functions [§ is similar to that of 
structure functions in DIS and leads to the introduction of generalized multi-parton cor¬ 
relation functions for fragmentation that parametrize 1/Q^ corrections, the moments of 
these correlation functions remain non-local quantities and a separation of short and long 
distances is not straightforward. Phenomenological hadronization models are at variance 
with the light-cone expansion of Ref. p, as they typically lead to 1/Q power corrections. 
The question of 1/Q versus 1/Q^ power behaviour is experimentally unsettled |^. 

The second moment of single-particle inclusive cross sections summed over all hadrons 
is of particular interest. Because of the energy conservation sum rule 


/ dxxDf(x,fj,) = 1, 


H 


the fragmentation functions disappear from the integrals 


^ 1 D dap ^ 1 D j 

ap = y - dxx —— = > - dxxCp\ 
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which can therefore be calculated in perturbation theory up to corrections suppressed by 
some power of A/Q. The normalization is such that the total cross section atot = o'l + ot- 
For the longitudinal cross section 


o'l — O'o 


— + (14.583 

71 


1.028Nf) 
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where = cngiQ) and Nf is the number of active fermion flavours. While dispersion 
relations relate the total cross section to the operator product expansion of a current 
correlation function, the power corrections to the longitudinal or transverse cross section 
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cannot be inferred from snch methods. From the theoretical point of view, the longitn- 
dinal and transverse cross sections can be considered as event shapes, and little has been 
known abont their power corrections nntil recently |^-|^. In Ref. [^] a 1/Q power 
correction to the longitndinal cross section was snggested as a conseqnence of phase-space 
redaction in the one-glnon emission diagram, when calcnlated with a massive glnon. At 
hrst sight, this conclnsion seems to be again in conflict with the expectation |P] that 
fragmentation fnnctions need to be corrected only at order 1/Q^. 

The present paper is devoted to a theoretical analysis of power corrections in frag¬ 
mentation processes.0 We collect resnlts of phenomenological character, which we hope 
can provide nsefnl gnidance in analysing fragmentation data collected at collid¬ 

ers and in extracting the strong conpling from the longitndinal cross section as well as, 
possibly, other event shape observables. The interpretation of onr resnlts in operator 
language will be given in a subsequent paper. 

In Sect. 2 we summarize the method, based on the infrared sensitivity of Feynman 
graphs, that allows us to estimate the x dependence of power corrections and clarify the 
assumptions that enter this approach. In Sect. 3, we give results for the power correc¬ 
tions to the quark and gluon coefficient functions. As will be explained, these power 
corrections factorize, so that the hnal estimate of power corrections to the fragmentation 
cross sections is given as a convolution with the leading twist fragmentation functions. 
A large part of Sect. 3 explores numerically the differences in various implementations 
of the method. Their comparison suggests an effective parametrization of 1/Q^ power 
corrections, presented in Sect. 3.5, which captures the gross features of the x depen¬ 
dence. This parametrization could be applied to LEP data iB§. as a substitute for 
phenomenological hadronization corrections obtained from Monte Carlo programs. We 
also hnd that the power expansion is singular in the region where the detected hadron 
is soft and emanates from a gluon jet. For small values of x, the effective expansion 
parameter is (A/(x(5))^. The consequences of these small-x singularities are pursued 
in Sect. 4, where we show that their resummation leads to a 1/Q power correction to 
the transverse and longitudinal cross section, thus resolving the apparent conflict with 
the light-cone expansion at hxed x. In Sect. 5 we argue that the large second-order 
correction in (^) is not accidental, but reflects the fact that the longitudinal cross section 
receives a 1/Q correction. This leads us to examine yet higher-order corrections in a 
certain approximation. We discuss their effect on a determination of a* and the energy 
dependence of the longitudinal cross section with an eye on generic features that would 
equally apply to other event shape variables. 

The results and conclusions of Sects. 3 and 4 overlap with those obtained inde¬ 
pendently by Dasgupta and Webber [0, who use a somewhat different model for the 


gluon contribution to the fragmentation functions. We discuss this difference and other 
ambiguities inherent in the method in Sect. 3. 


^ A preliminary account appeared in Ref. 0- 
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2 Method and assumptions 

2.1 Ultraviolet dominance 


Although the method we use to estimate the x dependence of higher-twist corrections 
has previously been used [^, ^ ^ in DIS, we find it useful to repeat the main ideas 
and to spell out the assumptions. Since the method applies without conceptual difference 
to fragmentation and DIS, it might be helpful to have DIS in mind as an example, for 
which the language of higher-twist corrections is familiar and higher-twist corrections 
can be interpreted within the operator product expansion. 

We start from the observation that the separation of leading-twist from higher-twist 
is not unique. This is not apparent in low-order perturbative calculations in the MS 
scheme. Imagine, however, that the factorization in transverse momenta that is im¬ 
plicit in were implemented by a rigid cut-off /x, such that only contributions from 
transverse momenta kt > ^ were included in the coefficient function Cp. Then one 
would find a term lnQ^//x^, whose cut-off dependence is cancelled by the /i dependence 
of leading-twist fragmentation functions and in addition power-like cut-off depen¬ 
dence, starting with which is cancelled by higher-twist contributions. Therefore 

the leading-twist contribution as a whole depends on the prescription used to implement 
the cut-off, just as the separation of leading-twist coefficient and fragmentation functions 
in d^) does. At first sight, such prescription dependence seems to be avoided in the MS 
scheme, because power-like dependence on the factorization scale /x does not exist. The 
problem reappears, however, because the coefficient function now has a factorially diver¬ 
gent series expansion in as (referred to as infrared renormalon divergence). Summing the 
series again requires a prescription and the prescription-dependence is power-suppressed 
precisely as the cut-off dependence above. In both cases, the ultraviolet renormaliza¬ 
tion of higher-twist operators must be performed consistently with the definition of the 
leading-twist coefficient function. The sum of leading-twist and higher-twist contribu¬ 
tions is then unique. For the case at hand, setting /x = Q in (|^), we can write the leading 
power ambiguity in the leading twist coefficient function as 


SC‘p(x) oc A‘ 2 ^p{x) 


( 6 ) 


times, possibly, logarithms of A/Q. The functions A 2 p{^) calculable in a certain 
approximation, as explained below. 

Whichever point of view one prefers, there are two immediate conclusions: first, from 
the phenomenological point of view, higher-order contributions in perturbation theory 
and higher-twist corrections are inseparable and should be described by one parameter; 
second, precisely for this reason, some information on higher-twist effects can be ob¬ 
tained from the infrared (or large-order) behaviour of perturbation theory. If we could 
enumerate the higher-twist operators as in DIS, this piece of information would refer 
only to the ultraviolet regularization properties of these higher-twist operators, with an 
X dependence given by A 2 p{x) in (|]). In the following we will assume that the x de- 
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pendence of the entire higher-twist contribution is proportional to A 2 p{x) and refer to 
this assumption as ‘ultraviolet dominance of higher-twist corrections’. Eq. (|^) is then 
replaced by 


dap 

dx 


[X 




dz 


CUz,QVf^")+K^2Ay{z) 




+KlA\p{z) ^ + • • • 


Di{x/z,ix). 


(7) 


The dots denote higher power corrections and we anticipated that only even powers in 
1/Q occur. The functions A\p{z) will be given in Sect. 3.5 and the are adjustable, 
^-independent constants that should be determined by comparison with experimental 
data. Note that here we slightly differ from previous formulations, where the constants 
are hxed in terms of an ‘effective coupling’ M, ^ or equated to a universal constant PO 


(hxed by the relation to the IR renormalon ambiguity). In our approach, these constants 
can depend on the factorization scale and also the order of perturbation theory to which 
the leading-twist coefficient function Cp has been calculated, since the added power 
corrections partly parametrize higher-order perturbative corrections as well. Naturally, 
the constants should be ‘of order 1’. We emphasize again that for i = q, a. sum over 
light quark flavours (u, d and s) is already included in Al^p{z). 

Arguments related to infrared cut-off behaviour or infrared renormalons have been 
used repeatedly [|T^] to determine the scaling of power corrections with 1/Q. In this case, 
no new information is obtained for DIS, where the power behaviour is known from the 
operator product expansion (OPE). The dependence of power corrections on xpj is not 
constrained by the OPE and is given in terms of multi-parton correlation functions, which 
are already too complex to be extracted from experiment. The ultraviolet dominance 
hypothesis provides tremendous simplihcations, as the unknowns are reduced to a few 
constants rather than functions. Results that followed from applying this hypothesis to 
DIS structure functions have been found to reproduce experimental results on the xpj- 
dependence of higher-twist contributions unexpectedly well |I9| , p0| , pT| . This empirical 
success encourages us to try the same idea for fragmentation processes. 

The assumption that higher-twist corrections are proportional to their cut-off depen¬ 
dence (or renormalon ambiguity) does not become correct in any limit and therefore the 
resulting x dependence should be considered as a model. One should not insist on this 
proportionality too literally. Rather, the idea is that if the ultraviolet contribution to the 
higher-twist correction varies rapidly with x, one may expect that the full higher-twist 
correction also does, up to some smooth function. This expectation is best motivated 
by an analogy. Suppose we calculated some quantity to order a” in perturbation theory. 
It is customary to vary the scale p to get an idea about the size of the next term in the 
expansion. Of course, the /i-dependence gets exactly compensated, but if it is large (or 
small) the /i-independent terms are also expected to be large (or small). In the same 
way, although renormalon ambiguities are ultimately cancelled and unphysical, their x 
dependence should be indicative of the x dependence of the full higher-twist contribution. 
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Related to the choice of scale is the fact that in certain kinematical regions the scale of 
the hard process turns out to be of the form e{x)Q, with e{x) typically given by a power 
of a; or of 1 — X, because of phase-space restrictions for gluon emission. The physical scale 
is thus parametrically smaller than Q, and the power expansion is naturally organized 
in terms of K/{e{x)Q) rather than A/Q. Since kinematic restrictions affect large orders 
in perturbation theory, the ultraviolet behaviour of higher-twist operators and the entire 
higher-twist correction equally, the method outlined here will reproduce such parametric 
enhancements. Such enhancements will be seen to be crucial in understanding behaviour 
of power corrections to integrated fragmentation cross sections. 

While the assumption of ultraviolet dominance of higher-twist corrections is dehni- 
tively crude and clearly misses some features of higher-twist corrections (such as op¬ 
erators that do not mix into leading-twist ones through their ultraviolet behaviour), 
it should perhaps be measured by the simplicity with which some generic features of 
higher-twist corrections (and higher-order corrections in perturbation theory) can be 
incorporated. 

2.2 Calculating A\^p{x) 

In this subsection we describe the calculation of the functions A\p[x) in ( 0 ). They can 
be obtained either from lowest-order Feynman diagrams with an infrared regulator or 
from infrared renormalons in the large-order behaviour of the coefficient functions Cp. 
The precise form of A\ p{x) depends on the method chosen, although some equivalences 
can be found ||22[. We adhere to the second method. In this case, the functions A\p{x) 
are dehned as the residues of IR renormalon poles of the Borel transform of Cp. This 
is still not practical, because it would require us to calculate the perturbative expansion 
of Cp to all orders. We therefore approximate the perturbative expansion by the series 
generated by inserting any number of fermion loops into the gluon line in diagrams that 
contribute at order 

This class of diagrams falls into subclasses according to which parton is registered 
in the hnal state and where it originates from. Fig. |l| shows a contribution to the 
partonic cross section da'p/dx, where the registered quark originated at the primary 
virtual photon vertex. Another ‘primary’ quark contribution is obtained if, in Fig. ^ a 
gluon line instead of a fermion loop is cut. The registered quark can also come from a cut 
fermion loop. We refer to this contribution to dap/dx, shown in Fig. ^ as ‘secondary’ 
quark contribution. Finally, the gluon cross section dap/dx is obtained from diagrams 
similar to that in Fig. but with a cut gluon line instead of a fermion loop. In addition, 
diagrams with counterterm insertions have to be included. The antiquark cross section 
is identical to the quark cross section. 

The set of diagrams with fermion loops is not really relevant by itself, as it typically 
yields small contributions to higher-order perturbative corrections. The main idea is 
to reinterpret the series in {NjasY given by this set of diagrams as a series in 
and to restore the full QCD /3-function coefficient /3o = —l/(47r)[ll — 2A^j/3] from the 
dependence on Nf ‘by hand’. This substitution seems difficult to justify, but comparison 
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with exact low-order results shows that it reproduces exact results quite well and that 
keeping corrections (/^oc^s)"' in higher orders resums important contributions. Further 
motivation can be found in |^, ^ ^ . 

The restoration of /do for the secondary quark contribution is especially delicate, since 
the fragmentation cross section depends on the ‘internal structure’ of the fermion loop, 
that is the phase-space over ki and /c2 (see Fig. |]) is non-trivially weighted. In such 
a situation, the association of fermion loops with running of the coupling is lost. To a 
certain extent, this instrinsic limitation of the method can be investigated by considering 
a fictitious theory with scalars rather than fermions in the fundamental representation. 
In this case we would compute scalar loop insertions and restore the full QCD /3-function 
coefficient /3o = — l/(47r)[ll —iVs/6] from the dependence on Ng. We would like to obtain 
the same result as with fermions, since in both cases the important contributions are 
related to the non-abelian contribution to /3o- For observables that resolve the internal 
structure of fermion or scalar loops, however, the resulting functions A\p{x) are not 
identical. This difference leads to some model dependence, which is analysed in Sect. 3. 

The fermion-loop approximation (and the massive gluon calculation, discussed in 
Sect. 2.3) also neglects multiple gluon emission. This is an obvious shortcoming of the 
method relevant at small values of x. The resummation of Inx contributions is crucial for 
particle multiplicities and the shape of fragmentation functions in the small-x region. In 
the following we consider the expansion in at fixed x, leaving aside the question of the 
small-a: asymptotics at fixed Q^. The order of limits is probably important. Renormalon 
and small-x resummations select different sets of diagrams and it is not known how to 
combine the two in a systematic way. See Sect. 4 for a further discussion of the problem. 
Related questions are also being discussed for Drell-Yan production [^, |2^ . 


The evaluation of the two sets of diagrams shown in Fig. and Fig. for an arbitrary 
number of fermion-loop insertions, and for their sum, is straightforward by means of the 
dispersion technique discussed in Refs. p3|, |^. For a moment, let us consider only 


diagrams with a cut fermion bubble. It is convenient to organize the calculation such 
that the integral over the gluon virtuality (invariant mass of the qq pair) is done last. 
The main object of our interest will be the distribution in the gluon virtuality To 
dehne it, consider for example the contribution of the diagrams in Fig. 1 or 2 to dap/dx, 
summed over the number of fermion loops. It is given by^ 


da 




dk"^ 


dx 



dap’^^ 

(x 

|i + n(F)|' 

dx 

(’ QV 


( 8 ) 


where n(/c^) is the MS-renormalized fermion bubble, and the mass distribution can be 
written as 

^ y'dLips[pi,p2,fc] {2'Kf5^^\q-pi-p2-k)^M^,yMl,^, 

^The following equation should be considered as schematic, because the limit ^ 0 needs some care. 
The correct treatment of this limit leads to (0 below. Alternatively, the equation can be understood 
as an expansion in Us- 
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Figure 1: Example of a ‘primary’ quark contribution to the squared amplitude for 
dap/dx. The set of all diagrams includes all attachments of the chain of fermion loops 
to the external quark line, the diagrams with a cut gluon line and an arbitrary number 
of fermion-loop insertions. 



Figure 2: Example of a ‘secondary’ quark contribution to the squared amplitude for 
dap/dx. The set of all diagrams includes all attachments of the chains of fermion loops 
to the external quark line and an arbitrary number of fermion-loop insertions. 




■ ^ fdUps[ki, fca] {2Tiy5^'^\k - ki - fcs) ^Tr[7^ ^17^/ ^2] ^P;’i'(2:) , (9) 

where ^ = k‘^/Q‘^. In Eq. (^, the notation for momenta and Lorentz indices corre¬ 
sponds to Fig. 2, and /3 q = Nf/iQn). We denote by [p] ([s]) the ‘primary’ (‘secondary’) 
quark contribution, while /(iLips[...] are Lorentz-invariant phase-space integrals, and 
is the matrix element for the primary 7* —> qqg amplitude squared, divided 
by the square of the quark electric charge. The projections are such that (in the 7* rest 
frame) 


V 


L-\-T 


{x) = 


Qiiii' 


X 


-pL^J 


X) = 


4|fci|2 


(5 X 


2kiq\ 

2 ) ’ 

(10) 

q J 


2kiq\ 
g2 ] 

(11) 


for the secondary contribution to the total and longitudinal fragmentation function. For 
the longitudinal projection we used q^Aif^u = 0, and the primary quark contribution 
is obtained by replacing ki with pi. The normalization factor 2/(31 ® may seem 

peculiar. It is chosen such that for inclusive quantities the distribution function coincides 
with the result that would be obtained from computing Os-corrections with a massive 
gluon. 

Note that in the case of the primary quark contribution the phase-space integral over 
ki, ^2 is proportional to so that the result takes the form of the one-loop diagram 
calculated with a gluon of mass k^. This ensures equivalence with the massive gluon 
calculation for this contribution. For the secondary quark contribution the projector 
depends on ki and modihes the phase-space integral. This is probably the simplest 
example of how event shapes in general are sensitive to the internal structure of fermion 
loops. 

Given the invariant mass distributions in the contributions of fermion-loop dia¬ 
grams to higher-order perturbative corrections are obtained in terms of the logarithmic 
moment integrals 


Jn{A^ [ 

JO 


in"e 


d dcTp’^^ 




( 12 ) 


d^ dx 

The sum of the series, dehned by a principal value prescription for the Borel integral, 
equals 


da 


(NNA) 


dx 


Jo df dx 


dap . , dap 


dx 


dx 


X, 0) 


(13) 


d^ dx 

where < 0 is the position of the Landau pole in the strong coupling and the function 
<h(^) is specihed in Eq. (2.25) of the second reference in |2^. We shall make use of these 
results in Sect. 5. 
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To obtain the functions Al^p(x), it is not necessary to perform the hnal integration 
over (^). The infrared renormalon residues can be read off as coefficients of non- 
analytic terms in the expansion of the /c^ distribution dap/dx{x,^) at small ^ 


dx dx 


0) + ^^ + ^2,p{x) ^In^ + Ai^p{x) 


e'lne + ... 


(14) 

Eq. (p!4D thus supplies all information required for the convolutions in (J^). 

Up to this point, we considered the partonic fragmentation cross sections H = q,g in 
(2). To extract the coefficient functions Cp, the partonic fragmentation functions Df^ 
have to be calculated and subtracted in the same approximation. The need for subtrac¬ 
tions is also reflected in the fact that the limit ^ 0 in ( [T^) does not exist, because 

dap/dx 


In^ at small As a consequence, the integrals in (|T^ and (|^ also diverge at 
small Once the partonic fragmentation functions are computed, the subtractions can 
be implemented into these equations by generalizing the case of ultraviolet subtractions 
discussed in to infrared subtractions. 

We emphasize that the power corrections added in (^) (that is the functions A„ p(a;)) 
depend on the factorization scheme and, for instance in the ‘annihilation scheme’, would 
differ from those in the MS-scheme, because the dehnition of the leading-twist fragmen¬ 
tation function can include an arbitrary set of power corrections. We will be working in 
the MS-scheme. In this particular scheme, we do not need to perform the subtractions 
explicitly and the functions An^p{x) can already be obtained from the partonic cross 
sections. To see this, let us sketch how collinear factorization is performed for the set of 
diagrams considered above. 

To solve d^) for Cp it is convenient to count powers oi Nf for a given diagram, where 
Nf denotes the number of flavours. We expand the partonic fragmentation functions 


D 


[X 


= DT 


[x 






(x) + ... 


(15) 


where Dx^"''{x) are power series in NfUg. The leading order (in 1/A^/) contributions are 

(16) 

(17) 


DI’A\ 

(x) 

= ^(1 - x) , 


(x) 

1 

1 

1 

£) 9 ,[ 0 ]| 

[x) 

= 0, 


[x) 

- "" Pn^nix] 

9 

27ie ^ ^ 

Nf/2 [x^ + 

(1- 

- xY] denotes 


1 - ’ 


x) + 0(Ny,) , 


(18) 

(19) 


wllGrG P ^ ^ -i_ f 1 TIto A P QT^l i-hhi-ri rr fn n r*-hi— 

Nf/(67r) is the fermionic contribution to the (3 function, and e = (4 — d)/2. Note 
that contains an entire series in NfOg, of which we show the hrst term only. 
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The factor 1/{1 — (ilag/e) in comes from counterterm insertions in fermion 

loops. Since the Born diagram leads to a non-vanishing coefficient function only for 
C'l’(x) = 5(1—x), we hnd for the coefficient functions, in the fermion-loop approximation, 


Cl{x) 



d(^L 

dx 


CpUs 4(1 — x) 
271 X 


( 20 ) 


Cl(x) Dfi , 

C-I(x) = 0 - _ DfA 


( 21 ) 

( 22 ) 


C^x) 


_ T^q,[i] _ r<9 ^ ng,| 

dx 


(23) 


The asterisk denotes convolution. The coefficient functions have hnite limits as e —>■ 0. 
Note that in the fermion-loop approximation, the longitudinal gluon coefficient function 
has no corrections beyond hrst order. The graphs with uncut fermion loops are cancelled 
by analogous contributions to H®. On the other hand, the longitudinal quark coefficient 
function and the transverse coefficient functions involve non-trivial subtractions. 

However, as long as we are interested only in power corrections (infrared renormalons) 
and not in the complete higher-order perturbative corrections, we do not need to cal¬ 
culate these subtractions explicitly. The partonic fragmentation functions needed in 
our approximation have convergent series expansions and do not contribute to infrared 
renormalon residues. This property is specihc to minimal subtraction schemes, where 
partonic fragmentation functions are expressed as pure poles in e. For the ^ distributions 
introduced above, this implies that the In,^ terms are removed, but all other non-analytic 
terms, and in particular the functions Hp()r), remain unaltered by factorizing the frag¬ 
mentation functions. 

It is important that the mixing between quark and gluon operators at leading twist 
repeats itself at the level of power corrections and renormalons. Although the gluon 
coefficient functions are convergent series in the fermion-loop approximation, this does 
not mean that gluon fragmentation does not give rise to power corrections. Because 
of mixing, the infrared renormalons in quark coefficient functions due to the secondary 
quark contribution are cancelled by ultraviolet renormalons in quark matrix elements 
of gluon operators of higher twist. The same gluon operators also have non-vanishing 
gluon matrix elements, which are set to zero only in the formal large-limit. 


2.3 Massive gluon scheme 

It can be shown that for sufficiently inclusive observables, the functions A\ p{x) 

obtained from the residues of IR renormalons in the above approximation to the per¬ 
turbative series can also be obtained from the corrections alone, provided they are 
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computed with non-vanishing gluon mass. In this case, these functions can be read off 
from non-analytic terms in the small-mass expansion.^ This equivalence does not hold 
for general event shapes, which are sensitive to the internal structure of fermion loops, as 
is the case for the secondary quark contribution in Fig. (and as discussed in Ref. ||14|| ). 

Although in this situation a formal justihcation of the massive gluon scheme is lacking, 
except, perhaps, as an ad hoc implementation of the general idea of scale separation by an 
explicit infrared cut-off, we hnd it useful to check whether the massive gluon calculation 
reproduces the gross features of the fermion-loop calculation. If so, one could take 
advantage of the relative simplicity of the massive gluon calculation. 

With this motivation in mind, we will present also the distribution functions in ^ in 
the massive gluon scheme (see also |^). We will hnd that the comparison between the 
fermion loop and massive gluon calculation has to be done case by case, with varying 
conclusions. 


2.4 Summary of distribution functions 

In this subsection we present our results for the distributions in ^ corresponding to the 
detection of the primary and the secondary quark, as well as those corresponding to 
the detection of a gluon, calculated in the massive gluon scheme. We give the exact 
expressions, as well as the small-.^ expansions corresponding to (|T^, and we extract the 
coefficients Aj p{x) of the power corrections. As explained above, the gluon coefficient 
function is trivial in the large-order approximation we have adopted and does not indicate 
power corrections. The role of the gluon coefficient is taken by the secondary quark 
contribution. Fig. Q, to the quark coefficient function. The results presented below refer 
to the sum of quark and antiquark contributions, which amounts to multiplying the 
quark contribution by two. 


2.4.1 Primary quark contribution 


For the primary quark contribution to dap/dx (Fig. [^), the distribution in coincides 
with the partonic quark fragmentation cross section calculated with a massive gluon. 
Dehning ^ = k'^/Q'^ as before we have 


1 dal^^ 
(To dx 


Cpc^s 

2ti 


■ 2-0(i-e 
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1 — X 


^)+|(2a; + 301n^ 


(1 




(1 - x){x + i) 


(24) 


1 

(To dx 


CpOis 

271 


■2- <^^(1-0:) 


2(1 + 0 " (Li2(-0 - i Iny + lnUn(l + 0 + y 


"^The restriction to non-analytic terms is clear. Analytic terms in also arise from large k^, where 
the propagator l/(fc^ -I- A^) can be Taylor-expanded in A^. 
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y 1 — X 
The expansion in ^ gives 
1 Cpa 


l_,_2^ + ^^+5Ch.> 


2(1-x) {1-xy 

+ 3^ 


X x'^ J (1 — a;)(a; + .^) 


(To dx 2n 


■ 2 ■ 


1 + ^ In^ Al’jl\x) + d - - 2 - , " + - In —^-- 

’ \x [1 — a;J+ X x(l — x) ^ 


( 25 ) 




(26) 


where 



= 2(5(1 - x) + -, 

X 

(27) 

Atl!{x) 

= (5'(1 - x) + 4(5(1 - x) + 

(28) 


The 6 functions arise from expanding singular functions around the phase-space bound¬ 
ary 1 — a; — For the transverse cross section, we keep only non-analytic terms in the 
expansion; then 


1 d(T^ 
tq dx 


CpOis 

271 


- In^ 


/ 1 + x"^ 

v[i-2^]+ 


+ 5^(1 





+ InO 

where 


AiP{x) = - 

4 

4 

+ 2 - 2(5(1-x) 

X 

+ 6'{l-x 

1 

' 
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dh 

+ 

At!f\x) = - 

2 


- x) 

1 

' 

1 

,H , 
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(29) 

(30) 

(31) 


2.4.2 Secondary quark contribution 

Because of cancellations between the longitudinal and tranverse cross sections it is more 
convenient to quote the distribution in ^ for the longitudinal contribution and the sum 
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of transverse and longitndinal contribntions, 
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-3«V)r(v^,i) , 


(To dx 




(^4^^ — a;^ + 2^x‘^ — ^"^x^ + a;^ + + 


6(1-0 


(a;^ + 3^) In^ 


where 


7 ^ (-20 - 2^a; - 20a: + x^ + A^x^ + 0a:0 In — 

C rv^Kj \ / /y»Zt 


24(l + a; + 0 


In a; In I - ^(1+0^ (a:^ + 3^) T(y^,a;) 


/•i , In(At) i r ^ , 1 — iA , A, 
= dt -—^ = — In A In-^— In — In 


T(A, x) = [ dt 

Jl/x 


1 /x 1 + A2t2 2A 


1 — lA , A , a; — lA 

--— In — In-- 

1 + iA X a; + iA 


+ Li2(iA) - Li2(-iA) - ^Li2 - Li2 . (34) 

The ^ 0 limits and the hrst non-analytic terms in the expansion at small ^ are given 
by 


1 (i(T|^’^^' CpCts 

(To dx 271 
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(36) 

(37) 

(38) 

(39) 

(40) 


Note that the expansions in ^ are valid only for x > a/^, although x can be as small as 
The minimal invariant mass of the qq pair, however, is attained at x = a/^. 

The ^ > 0 limits in (|35D and (^) require explanation. First note that although the 

diagrams in Fig. and Fig. start at order a'^, our resulting distribution functions are 
written as order as- Technically, this can be best understood by considering the Borel 
transform of the series generated by these diagrams. Since there is no diagram at order 
as, one would expect the Borel transform to vanish linearly in u, the Borel parameter, at 
small u. However, the g ^ qq splitting amplitude contains a collinear divergence, which 
manifests itself as a pole 1/m. Consequently, the Borel transform approaches a constant 
for small u, which corresponds to an ag contribution. The collinear pole in u indicates 
that collinear subtractions need to be performed. Indeed, if Cl * is subtracted as 


required by (|2T|) , the constant term in the Borel transform vanishes and the secondary 
quark contribution to the quark coefficient function is of order a^, as it should be. 


To check this cancellation, we observe that the leading term in (^) can be rewritten 


as 


CpOis 

2tt 


■ 2 ■ 


— 6x + 2x^ + 6 Inx 


X 


= 2 ■ 


3 

ITj 


■ [Cl * p. 


9^q\ 


X 


(41) 


with Cl dehned 


m 


As expected from the above discussion, the ^ 0 limit of 


the secondary quark contribution is just the gluon coefficient function convoluted with 
collinear splitting into a qq pair. The factor 2 accounts for the sum over quark and 
antiquarks. The factor 3/Nf is the Isnge-Nf limit of as/{2'n') ■ InQ^/A^ associated with 
the g ^ qq amplitude. 
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The In ^ term in the sum of the longitudinal and transverse cross sections can likewise 
be interpreted as a convolution of two subsequent q{q) —>• g and g —>• q{q) splittings, as 
indicated in (pbl) . 

We should stress again that these complications associated with collinear factoriza¬ 
tion do not affect the non-analytic terms in the expansion in except for In^, and 
therefore do not affect the expressions for p. Note that for heavy quark fragmenta¬ 
tion, the secondary quark contribution would not give rise to power corrections, because 
the invariant mass distribution is cut off at drug. 


2.4.3 Massive gluon scheme 


For completeness and later comparison, we quote the distribution functions in the mas¬ 
sive gluon scheme. The distribution functions are given as the partonic fragmentation 
cross sections to order ag, computed with a massive gluon. The quark fragmentation 
cross sections are identical to our primary quark contributions. Sect. 3.1.1. The gluon 
cross sections read 
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(To dx 
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l + ^-x 
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(42) 
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X' ' 4 ^ 

As ^ —^0, the longitudinal cross section reduces to (^0[) . Expansion in ^ results in 
1 daf^^ 


(43) 
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where again we have kept only the leading term and the first few non-analytic terms. 
The coefficients of power corrections are now given by 
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,MG, 
2,L ' 
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( 47 ) 




a:-’ 


X'^ 


^ 2 ,l+t{.x) = -+ 4 - 25(1 - x) , 

’ X 


(48) 


^i’X+ri'^) - - 85(1 - x) - 5'(1 - x) 

X 


(49) 


Our results in the massive gluon scheme coincide with those obtained by Dasgupta and 
Webber 


3 Power corrections to fragmentation functions 

In this section we analyse the x dependence of power corrections. We present predictions 
in the fermion-loop approximation and then discuss in detail the ambiguities related 
to the restoration of the gluon fragmentation component. To this end, we compare 
the fermion-loop calculation with the massive gluon calculation and check the effect 
of replacing fermions by scalars. Once again, ‘massive gluon calculation’ refers to the 
identification of power corrections through non-analytic contributions in the expansion 
of order-a^ corrections in a small gluon mass, as done in Ref. |^. Finally, in Sect. 3.5, 
we abstract from the discussion the generic dependences on x and formulate a simple 
phenomenological parametrization of 1/Q^ corrections consistent with our results. 


3.1 Results 

To illustrate the magnitude of the leading l/Q"^ power corrections to the fragmentation 
cross sections in e’''e“ annihilation in the fermion-loop approximation, we rewrite m as 


dap 

dx 


ix,Q^) = Fp{x,Q‘^) 


where 


1 -F H2,p{x, — + . 


(50) 


Fp{f Q") =T.\Cp* AI (:r, Q^)=Y.I - CUf QVf") Di{x/z, /i) ( 51 ) 

denotes the leading-twist cross section. The power correction is given by 

p * Di] (x, Q'^) 


Fp{x,Q^) 


(52) 

In the following, we set iF| = 1- Since in the fermion-loop approximation we only have 
a quark (antiquark) contribution, no information is lost, as the over-all scale is then set 
by in (^). We use the ALEPH parametrization |P] of the (light) quark (and later. 
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gluon) fragmentation function at Q = a/s = 22 GeV and their value as(22GeV) = 0.164. 
In evaluating Fp{x,Q‘^), we use the lowest-order approximation to the leading-twist 
coefficient function Cp. 

The result for H 2 ^p{x) for the longitudinal and transverse fragmentation cross sec¬ 
tions, summed over all hadrons, as a function of energy fraction x is shown in Figs. ^ and 
respectively. The sum of longitudinal and transverse cross sections, which differs from 
the total cross section by the small asymmetric contribution due to interference, 

is shown in Fig. We recall that Fp/Fp ~ so that the relative magnitude of power 
corrections is enhanced for longitudinal fragmentation as reflected by the scales on the 
vertical axes in Fig. ^ and |[ We note that for x > 0.2 the power correction to both, 
longitudinal and transverse, cross sections is dominated by primary quark fragmentation. 
Such a qualitative behaviour is expected, since the average energy of a quark connected 
to the ‘hard’ 7* or vertex is much larger than the average energy of a quark orig¬ 
inating from gluon splitting. Since for the primary quark contribution the restoration 
of Po from the dependence on Nf is unproblematic and the fermion-loop and massive 
gluon calculations coincide, the method yields an unambiguous prediction at x > 0.2. 
At smaller values of x the secondary quark contribution becomes important for the lon¬ 
gitudinal cross section and for the sum of longitudinal and transverse cross sections, 
but remains small for the transverse one. Gomparing Figs. with the corresponding 
figures in Ref. |^, obtained from the massive gluon calculation, one concludes that the 
results are qualitatively similar for the longitudinal cross section, but differ drastically 
at small x for the others. The difference can be traced to the absence of a 1/x^ term in 
Eq. (^81), as compared to Eq- (HD, and will be discussed in detail 

in Sect. 3.4. 

Let us comment on the normalization of H 2 ^p{x). With = 1, the overall normal¬ 
ization of A 2 ^p was obtained by setting CpUs/{27i) ^ In ^ —> IsP/Q"^ in the distribution 
functions. If we took the normalization literally from the IR renormalon ambiguity of 
the Borel integral, divided by tt, we would have to substitute 


Cpas 

27r 


,, , Cpe^/^ A2 ^^A2 
^ ^ 27r(-/?o) ~ 


(53) 


where A = Q exp(l/(2/3oas(g))). With this substitution, using as as given above and 
H 2 . l+t{x) ~ 15 for intermediate x, the power correction in Fig. ^scales as (0.7 GeV/Q)^, 
very similar in magnitude to what is obtained with an ‘effective coupling’ [|T^. For 
phenomenological analyses we suggest that the overall normalization be fitted to the 
data. 


3.2 Restoring gluons 

In the above analysis, all power corrections are obtained as convolution with the quark 
fragmentation function. But the original motivation was to compute fermion loops only 
to trace the dominating non-abelian (gluon) contribution through the dependence on 
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Figure 3: Shape of 1/Q^ power correction H 2 ^l{.x) to the longitudinal fragmentation 
cross section. Dashed line: primary quark contribution. Dotted Line: secondary quark 
contribution. Solid line: sum of both. 
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Figure 4: Shape of 1/Q^ power correction H 2 ^t{x) to the transverse fragmentation 
cross section. Dashed line: primary quark contribution. Dotted Line: secondary quark 
contribution. Solid line: sum of both. 
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Figure 5: Shape of 1/Q^ power correction H 2 ^l+t{x) to the sum of longitudinal and 
transverse fragmentation cross sections. Dashed line: primary quark contribution. Dot¬ 
ted Line: secondary quark contribution. Solid line: sum of both. 

Nf. Eventually, the secondary quark contribution should therefore be reinterpreted as a 
gluon fragmentation contribution. 

Such a reinterpretation faces ambiguities, and we would like to single out the gross 
features, which can be considered as unique. To have a basis for comparison, we hrst 
‘deconvolute’ the g ^ qq splitting from the secondary quark contribution. We dehne the 
‘effective gluon’ coefficient functions Af through 

* Pg^q] (a^) = (54) 

which should be convoluted with the gluon fragmentation function, replacing A^p con¬ 
voluted with the quark fragmentation function. The superscript g ^ q indicates that the 
present effective gluon coefficient has been obtained by deconvoluting a quark emission 
amplitude. Below we will consider an analogous coefficient obtained from a scalar quark 
emission amplitude. Although 

* D, = * D, (55) 

holds only to leading logarithmic accuracy in Q^, this is the closest one can get to 
reinterpreting the secondary quark contribution as a gluon contribution. The functions 
yyi x) can be compared directly with the Af’p*^(a;) obtained in the massive gluon 
calculation. 

To test the sensitivity of the deconvolution to the particular structure of the quark- 
gluon vertex, we compute p in a (hctitious) theory with scalar particles rather than 
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fermions (quarks). Then we deconvolute the gluon-scalar splitting function Pg^sq{x) = 
a:(l — x)/2 to obtain the effective gluon coefficient The gluon interpretation 

of the secondary quark contribution is justifed only if this function, after convolution 
with the gluon fragmentation function, leads to results compatible with those obtained 
starting from quarks. 

To obtain the functions for scalars, one replaces (neglecting terms that vanish 
when contracted with the trace 


2 N f . , 

^ ■ -if Trbi. hlw h] 
Po ^ 




(56) 


in i) by 


2 N 

^ ■ “ir (^1 - k2)u{ki - k2)u' = 247r ■ Akiukiy^. (57) 

bo ^ 

When the phase-space over ki and k 2 is integrated unweighted, the integrals of the previ¬ 
ous two equations coincide. Because of the projection Vp-^ggi in (||), however, Eqs. (^71)- 
(^) for the secondary quark contribution are replaced by 



^2,L+r 


(x) 


5X3 g. 


- - + 2 - -x^ 


+ 


35 x 3 ' 53.3 5 + 35^ ’ 


5x3 2 ; 


- - + 2 - -x^ 


(58) 

(59) 

(60) 


^4,L+r 


(x) 


9 4 1 2 2, 

3^ 


(61) 


for scalar particles. Note that = 0. 

For the interesting case of a gluon splitting into two gluons, the expression replacing 
(^) is still a linear combination of only two Lorentz structures ki^ykiy and k^Qui,'. There¬ 
fore all ambiguities related to restoring are already exhausted by a linear combination 
of the results obtained for fermions and for scalars. In particular, if one computes the 
graphs with one cut gluon loop in an axial gauge, the coefficient of ki^k^,/ is independent 
of the gauge-hxing vector. 


3.3 Discussion: longitudinal fragmentation 


We now consider the secondary quark contribution to longitudinal fragmentation in 
more detail. First, we compare the secondary quark contribution (see Fig. ^ with the 
corresponding contribution in the massive gluon scheme. The ratio 


R{x) 


At^ * py) 

Air * D,(x) 


(62) 
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Figure 6: Ratio of secondary quark contribution in the fermion-loop approximation and 
gluon contribution in the massive gluon scheme. 



Figure 7: alix) (solid), ai^x) (dashed) as dehned in the text, compared with the 
massive gluon scheme (dotted). 


22 






is plotted in Fig. |^. Note that the fermion-loop approximation and massive gluon scheme 
coincide for the primary quark contribution, while the secondary quark contribution is 
smaller by a factor of about 3 than its massive gluon analogue. The difference in absolute 
normalization would be meaningful only if we tried to £x the constants in terms of a 
single (universal) number as in [^. If the are considered as adjustable parameters 
as proposed here, the important message conveyed by Fig. ^ is that the ratio R is flat, so 
that the same x dependence is obtained with the two methods. The remaining variation 
is insignificant within the experimental uncertainties in the parametrization of leading- 
twist fragmentation functions and the theoretical uncertainties of our method. 

Second, we compare the effective gluon coefficients obtained from secondary quarks 
and scalars as described above. The coefficients are given by® 


.n^n, N 72 lira: 


X 


15 39 ^ r^i ■ 

--V^^cos -lira: 

2a: 14^^ L 2 


3^7 ^ . ry7, ■ 

-V3:sm -In a: 

2 1-2 


AlTix) = 


8(1 - a:) 


x-^ 


{2 + 2x 


X 


(63) 


(64) 


for quarks and scalars, respectively. We recall that the gluon coefficient in the massive 
gluon scheme is simply 8(1 — x)/x^, see (^^. To compare the x dependence on a 
reasonable scale, we dehne 


A’Cl-M = ■ I ■ alix) , (65) 

’ x-^ 8 

3lfl‘“(^) = 54^ . 1 . . (66) 

The normalizing factors are chosen such that ai;,(0) = 1. The x dependence of is 
shown in Fig. ^ for all three procedures (massive gluon, fermion or scalar loops). We 
observe again very little x dependence of the residual functions aii^x) and conclude that 
the method yields a unique result for l/Q^ power corrections to longitudinal fragmenta¬ 
tion. The residual different x dependence is dehnitely beyond any accuracy that can be 
expected from the model. 


3.4 Discussion: transverse fragmentation 

For the transverse fragmentation function, the effective gluon coefficients read 

In a: 


Aiy{x) = -18- 


X 


51 15 ^ rv^, 

— +12-—Vascos — Ina: 
2a: 2 l 2 


®To compare absolute normalizations, the effective gluon coefficient obtained from quarks has to be 
divided by 3 and that obtained from scalars by 12 to take care of additional factors as in (il)- 
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27 r ■ \ ■ 

—- mx 

L 2 




= 0 . 


(67) 

( 68 ) 


We have written the result for the transverse part separately, rather than the sum of 
transverse and longitudinal (referred to as ‘total’ in this subsection), in order to empha¬ 
size the absence of a 1/a:^ term. Being absent for quarks and scalars, it follows that it is 
absent for gluon loops as well in any gauge; therefore, the statement is independent of 
how precisely the non-abelian contribution is reinstated. On the other hand, cancellation 
of the leading 1/x^ terms is most likely specihc for the approximation of considering one 
chain of fermion bubbles and, in this sense, accidental as a pecularity of the box graph. 
For 1/Q^ corrections, the leading term for small x, which is 1/x^ in this case, is indeed 
present for both transverse and longitudinal coefficients, see in Sect. 2.4.2. Since 
in large orders of perturbation theory diagrams with two chains of fermion loops are 
not suppressed, the cancellation that we observe does not imply, strictly speaking, any 
parametric suppression of the leading power correction for transverse fragmentation. A 
numerical suppression can hold, however, and it would be interesting to see whether it 
can be inferred from experimental data. Such a numerical suppression would be natural 


for the effective coupling approach of [12, 19 


Comparing these results with the massive gluon scheme, we observe that in this 
scheme the leading term l/x^ at small x is cancelled in the sum of longitudinal and 
transverse coefficients, see (^). This cancellation appears to be unphysical and related 
to the fact that the structure of the massive gluon fragmentation cross section (^) 
is over-simplihed. Indeed, the analytic terms in the ^-expansion of (|43|) all behave as 
^n/x'in+i^ On the other hand, non-analytic terms can arise only from the coefficient of the 
In l/(4^)-term, which does not have a non-trivial expansion in All terms In^ vanish 
for n > 3 in this scheme, another indication that the massive gluon scheme may not 
capture the generic x dependence of power corrections in the case of total fragmentation 
functions.^ 

We compare the x dependence of the effective gluon coefficients obtained from quark 
and scalar loops for the total fragmentation cross section in Fig. p. The relative normal¬ 
ization is chosen such that the asymptotic behaviours at small x coincide. Compared 
to Fig. 1^, we have not scaled out the small-x behaviour, because we want to emphasize 
that the qualitative differences at large x will be insignihcant for practical applications 
after convolution with the gluon fragmentation function, because primary quark frag¬ 
mentation dominates in this region. Opposite to longitudinal fragmentation, the large-x 
behaviour of tranverse secondary quark fragmentation is not unambiguously predicted 
by our model. While the scalar result vanishes linearly at x = 1, the quark result 
approaches a constant. 


®We discuss in Sect. 4 that the presence of power corrections for the total fragmentation 

function does not conflict with the absence of l/Q power corrections to the integrated gluon fragmen¬ 
tation cross section. 


24 












X 


Figure 8: Comparison of effective gluon coefficients obtained from quarks (solid) and 
scalars (dashed) for the sum of longitudinal and transverse fragmentation. 


Summarizing this discussion, we recommend the usage of the same functional form 
for the gluon coefficient in total fragmentation and in longitudinal fragmentation, but 
with a separate adjustable multiplicative constant (recall that for scalar quarks, in the 
present approximation, the two coefficients are in fact equal [Eqs. (|5^ ) and (^)], while for 
fermionic quarks they differ only slightly [Eqs. (l37|) and (|3^) ]). If the adjustable constant 
for total fragmentation turns out to be close to the one for longitudinal fragmentation, it 
would indicate that the cancellation of 1/x^ terms for transverse fragmentation observed 
in our approximation is operative to a degree. If it turned out much smaller than for 
longitudinal fragmentation, this would provide support for the cancellation observed in 
the massive gluon scheme. In addition, we may add a constant to account for the x —>■ 1 
behaviour. However, the value of this constant is not important in practice, unless it 
were abnormally large. 


3.5 A simple parametrization 

We summarize the analysis of this section in the form of a simple parametrization of 
power corrections, which we suggest to apply to the analysis of fragmentation data. We 
write the fragmentation cross section (^ as 

^(x,Q2) =-^(x,Q 2) + —--(x,Q^). (69) 

ax ax ax 
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The leading-twist cross section is as in (^. The 1/Q^ power corrections are parametrized 
as 
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(71) 


where Di denotes the leading-twist fragmentation function for parton i to decay into any 
given hadron, and the plus distribution is defined as usual, 


dx 


f{x) 


[1 


x\ 
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= J dx 
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f{x)-f{l) 


X 


(72) 


for any test function /. The coefficients in front of quark fragmentation functions are 
taken from Sect. 2.4.1. The coefficients in front of gluon fragmentation functions follow 
from the discussion in Sects. 3.3 and 3.4. Note that for longitudinal fragmentation, we 
have neglected the residual x dependence displayed in Fig. and chosen the simplest 
functional form. Let us add the following remarks: 

(a) The overall scale IGeV/Q {Q in GeV) has been chosen to set the overall mag¬ 
nitude for the constants Ck and d. With this scale, we expect these constants to be of 
order 1. 

(b) The constants Ck and d should be htted from the data. As mentioned above, 
the constant d is actually insignihcant, so that the parametrization depends on four 
constants only. It is important that these constants should be quoted only in connection 
with a value for /i chosen for the factorization scale in the leading-twist contribution and 
also together with the order of perturbation theory to which the leading-twist coefficients 
Cp have been used. In this sense, the constants should be considered /i-dependent. We 
recall that the fermion-loop approximation considered in this paper suggests Cg^i ~ Cg^p+T 
while in the gluon mass scheme Cg^p 3> Cg^ppp- 

(c) Since we did not consider additional logarithmic Q dependence of power correc¬ 
tions, the scale /i in Di above is undetermined. Without additional information, it is 
natural to set fi equal to the scale used in the leading-twist part of the cross section. 

(d) If the above parametrization of power corrections is used, no Monte Garlo cor¬ 
rection for hadronization should be applied in addition. 

(e) The ansatz can only be used as long as a; > A/Q, where A ~ 1 GeV is a typical 
QGD scale. At smaller values of x, higher power corrections such as 1/Q^ become as 
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important as l/Q^ corrections. In fact, the 1/Q^ expansion breaks down and would 
have to be resummed. In other words, (|70[) and (|7TD are corrections to the leading-twist 
prediction. If the power correction becomes as large as the leading-twist result, the 
ansatz fails. 


4 XjQ corrections to integrated cross sections 

For X > K/Q and not too close to 1, the power corrections discussed in the previous 
section arise only from regions of integration where two partons are collinear. In this 
section we consider moments of the x distributions, specihcally the second moment (^), 
related to the integrated longitudinal and transverse cross sections, which are infrared 
hnite. 

The moments can have qualitatively different power corrections from those at hnite 
X, because they include new infrared sensitive regions x —>■ 0 and x —>■ 1, related to soft 
partons. For small x this possibility is evident from the increasingly divergent terms 
l/{Q^x‘^)"' found in the secondary quark or gluon contributions. Thus, the integration 
over X must be performed before expansion in In other words, an inhnite series of 
power corrections in 1/Q'^'^ has to be resummed to obtain the integrated cross section. 
Specihcally, for the gluon contribution in the massive gluon scheme 

(73) 

for any n. Recalling that ^ should be interpreted as we see that the small-x 

region can produce a 1/Q power correction, parametrically larger than at hnite x. 

In Sect. 4.1 we collect the distributions in ^ for the various contributions to the 
integrated cross section. Sect. 4.2 contains some observations on the relation of small-x 
behaviour and 1/Q corrections. 


4.1 Longitudinal and transverse cross section 

The primary quark contribution to the longitudinal and transverse cross sections is given 
by 
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The primary quark contribution does not contain odd powers of and therefore no 
1/Q power correction. 

Because of large cancellations, it is again useful to present the distributions for the 
longitudinal and total secondary quark contribution. We have 
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where 
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and no further half-integer powers of ^ appear in ([76D beyond 0'^^- total secondary 
quark contribution is simply 
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In the massive gluon scheme, the gluon contribution, to be compared with the secondary 
quark contribution, reads 
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Both the secondary quark contribution and the gluon contribution in the massive gluon 
scheme contain a term, although with different coefficient. This term cancels in the 
sum of the longitudinal and transverse cross sections. Summing atot = cTi^^ + crT^^ + cTi^T 
we reproduce the distribution function for the total cross section given in Appendix B of 
Ref. Note that because the longitudinal and transverse cross sections are infrared 
finite, we continued to neglect explicit infrared factorization. However, the secondary 
quark contribution now contains the order-a^ gluon contribution. 


4.2 Discussion 

We now discuss some general aspects of how 1/Q corrections are generated in the gluon 
(secondary quark) contribution to the integrated cross section. As far as these general 
aspects are concerned, the fermion-loop approximation and the massive gluon scheme 
are the same and we choose the latter for the discussion, because the analytic expressions 
are simpler in this case. 

We noted above, in Eq. ([75|), that 1/Q corrections are made possible, because the 
expansion parameter at finite but small x turned out to rather than To obtain 

the coefficient of in the integrated cross section, we need to pick up the most singular 
term in 1/x at every order in S,- Once this is realized, the real question is not why 1/Q 
corrections arise in the integrated cross section, but why they do not arise in the sum 
of the longitudinal and transverse cross sections. Indeed, the most singular terms in the 


expansion of the total gluon cross section (^31) are 
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and the general structure of the expansion is the same as for the longitudinal and trans¬ 
verse cross sections separately. Integrating this expansion term by term, we obtain, for 
the coefficient Oi of in the expansion of cr£’.^, the expression 
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= lim- |1 - 2 fi(-l/ 2 , £, 1 + £; 1)1 = 2(1 - ln2), 

e^O e 

so that Oi = 0, as is evident from ([78|) . From this exercise we see that the cancellation of 
1/Q corrections in the total gluon contribution is a property of the power expansion at 
hxed X only when all orders are included. Had we truncated the sum over n, a non-zero 
result would have been obtained. Furthermore, it is necessary to consider also analytic 
terms in the ^ expansion at hxed x to obtain the cancellation. This means that while 
1/Q corrections are tied to soft partons, collinearity is not important. 

We would like to note that these properties are reminiscent of those that lead to 
a cancellation oi 1/Q corrections in the Drell-Yan cross section in the same (fermion- 
loop) approximation [Q. Contrary to the Drell-Yan process, however, we know for 
fragmentation that 1/Q corrections must cancel after summing over contributions from 
all partons, because the total annihilation cross section receives at most 1/Q^ power 
corrections (leaving aside issues related to parton-hadron duality). On the other hand, 

looks rather 


L+T 


+ a' 


L+T 


the cancellation of 1/Q terms in cr£ -|- crf> and 1/Q^ terms in a' 
accidental from the diagrammatic point of view. 

The fact that the power expansion at hxed but small x runs in reminds us that 
the scale A of non-perturbative ehects must be compared not with the total energy Q 
of the process, but with the energy {Qx)/2 of the detected hadron. To see this more 


transparently, consider the longitudinal gluon contribution (^^. We may simplify the 
factor (1 + 0(1 + 'C lo 1, because it does not alter 1/Q corrections, and we omit the 
factor CpOis/{/2 tt). Now x times (^^ depends on ^ and x only in the combination x/\f^. 
The Borel transform B[F] of the perturbative expansion generated by cut and virtual 
fermion-loop diagrams is related to the ^ distribution F by the Mellin transform 


B[F]{u,x) = - 5 !!^ / dU-Z/pi+x) 

Tiu J at 




(84) 


Because of the dependence on x/y/^ only, the x dependence factorizes in the Borel 
transform as 
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with an ^-independent function B{u), which is easily obtained from (|4^) and (l8^). Its 
precise form will not be needed, but we note that B{u) is analytic for |m| < 1 with a pole 
at M = 1 related to 1/Q^ power corrections. The (formal) Borel representation is now 
given by 
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The X dependence is completely accounted for by setting the scale of the coupling con¬ 
stant to Qx/2, the energy of detected particles in the hnal state. Moreover, if we adhere 
to a strictly perturbative approach, the Borel integral diverges at inhnity for x < 2A/Q 
due to the Landau pole in the (one-loop) running coupling. Sensitivity to m = cxo in¬ 
dicates that the entire power expansion in 1 /Q^ needs to be resummed and that the 
renormalon approach is inapplicable at such small x. 

Let us now dehne moments in x with a cut on x that eliminates the small-x region: 


a- 
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Eq. 


gives 




1 _ ™1+7-2m 
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( 88 ) 


1 -(- 7 — 2m 

Taking Erst Xc = 7 = 0, we see that the 1/Q correction to the integrated longitudinal 
gluon cross section (corresponding to a pole of the Borel transform at m = 1/2) again is 
an immediate consequence of the scale being Qx. If the scale were Qy/x, the denominator 
would read 1/(1 — u) instead of 1/(1 — 2m). 

Keeping 7 = 0, we note that the pole at m = 1/2 is eliminated no matter how small 
Xc, as long as it is non-zero. Since in any experimental situation a minimum energy cut 
is required, the presence of 1/Q corrections in measured quantities depends on the value 
of Xc. If Xc ~ A/Q, the Borel transform is sharply peaked at m = 1/2 and the structure 
of large-order perturbation theory is indistinguishable from a 1/Q IR renormalon for 
practical purposes. 

Taking Xc = 0, we note that the position of the pole depends on the moment. In 
particular, depending on 7 , it is not bound to be integer or half-integer, a situation that 
never occurs in deep-inelastic scattering. It is this dependence on the precise weighting 


of long-distance regions (small x), first noted in |]^, that makes 1/Q power corrections 
to a^ and event shapes in general difficult to interpret in terms of operators (as in DIS) 
in the conventional sense. It also prompts some caution regarding the hnal conclusion 
on 1/Q corrections. It is conceivable that resummation of logarithmic terms In” x in 
perturbation theory modihes the approach to the small-x region. If, for illustration, 
resummation would turn x“^“ into in (^), a fractional power-correction would 

accordingly be obtained from (j^). However, at present we do not know how to combine 
small-x resummation and constraints from angular ordering with a renormalon analysis. 

The resummation of leading In x can formally be achieved by a modihcation of the 
DGLAP evolution equation, which replaces the fragmentation functions Di{x/z, Q^) with 
Di{x/z,x‘^Q‘^) under the convolution integral m- The evolution equation then reads 




asiQ) 
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Trading as{Q) for l/(—/3olnQ^) and writing the fragmentation function formally as a 


31 










Mellin integral 


Dg{N,Q^) = / , 

0 


(90) 


one can get the solution 

D{N, u) oc + 1 - iV)2T'33/^(2w - , (91) 

where 'ygg = 2Nc/{N — 1), b = —Airf^o and G{N,u) is analytic at iV + 1 — 2m > 0. 
We approximated the gluon-gluon splitting function by its dominant term at small 
X. As N ^ 1, the singularities at w = 0 and u = {N — l)/2 merge, and produce 
Dg{N,u) oc (1/m) exp[—2iVc/(&w)], which gives rise to the correct small-x behaviour. On 
the other hand, (0) exhibits singularities at m = {N -1-1- k)/2, with A; = 1, 2,..., which 
might indicate non-perturbative effects through power-suppressed ambiguities. Such an 
interpretation is, however, speculative, and we will not pursue it further in this paper. In 
general, one may ask, as for Drell-Yan production [^, ^ 
of this type also give information on power corrections. 


whether evolution equations 
In the above example, keeping 
only leading-lnx terms in each order of perturbation theory gives 2B{0) as residue of the 


pole at M = 1/2, while the correct result is 25(1/2), see (| 88 |) . In general, one has to be 
aware that systematic procedures to sum Inx terms may not work to power accuracy. 


5 Perturbative corrections to ai and the determina¬ 
tion of as 

A measurement of the integrated longitudinal cross section could eventually yield a 
precise determination of the strong coupling. To approach this goal theoretically one 
has to control higher-order perturbative corrections, as well as non-perturbative effects, 
both of which are expected to be much larger for than for the total cross section atot- 
In this section we consider perturbative corrections to written as 
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TT 
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with j3o = —l/( 47 r)[ll — 2Nf/3] as before. We approximate the exact higher-order co¬ 
efficient by its value in the ‘large-/3o’ limit, where /3o is restored from the term with 
the largest power of V/ at each order. Given the ^ distributions for the primary and 
secondary quark contributions to ctl in Sect. 4.1, the coefficients dn are obtained numer¬ 
ically from the logarithmic moment integrals given in Eq. (|I^), as discussed in detail 
in [^. The ‘large-/do’ approximation, called ‘naive non-abelianization’ (NNA) in 
reduces to the BLM scale setting prescription 
rewrite the exact ai correction in (^ as 


28] for M = 1. To see how it works, we 


di = 6.17- 0.7573/(-/do). 


(93) 
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n 

(MS) 

(MS) 


<i».W (V) 

(V) 

0 

1 

1 

0.036 

1 

1 

1 

11/2 

13/2 

0.052 

29/6 

1.121 

2 

29.82 

45.97 

0.060 

27.07 

1.203 

3 

164.1 

369.0 

0.064 

188.8 

1.226 

4 

944.1 

3441 

0.066 

1634 

1.219 

5 

5829 

3.73410^ 

0.068 

1.703 10^ 

1.210 

6 

3.940 10^ 

4.682 10® 

0.070 

2.088 10® 

1.205 

7 

2.948 10® 

6.70710® 

0.071 

2.95410® 

1.202 

8 

2.44710® 

1.086 10® 

0.073 

4.75110^ 

1.201 


Table 1: Perturbative corrections to ai as obtained from ‘naive non-abelianization’. The 
fourth column shows successive values for (JL,n/(^tot at Q = Mz and with as{Mz) = 0.118. 
The last two columns show a comparison of NNA coefficients in the V scheme {C = 0) 
obtained from fermions or scalars. The asymptotic ratio of scalar and fermion coefficients 
is 6/5; atot/o'o = 1-04 has been used. 

With —/do = 0.61 for Nf = 5, neglecting the second term gives an accuracy of about 25%. 
In particular, the approximation predicted correctly the large size of the second-order 
correction. 

We have calculated the coefficients dn in higher orders, in the MS scheme. The 
‘primary’ and ‘secondary’ quark contributions, and d^^^\ respectively, add to dn as 
dn = + 2d^%l/3. Columns 2 to 4 of Tab. |I| show the primary and secondary quark 

coefficients, together with successive finite-order approximations to <JL/atot based on 
these coefficients. The perturbative coefficients grow rapidly, especially for the secondary 
quark contribution. The hxed-sign growth and faster growth for the secondary quark 
contribution are directly related to an IR renormalon, indicating a 1/Q correction to 
secondary quark fragmentation as discussed in Sect. 4. 

Since the longitudinal cross section is not fully inclusive with respect to gluon splitting 
into a quark-antiquark pair, restoration of the non-abelian piece is not unique, just as it 
was not unique for power corrections. The resulting ambiguity presents a major difficulty 
for the extension of ‘naive non-abelianization’ (and BLM scale setting) to hadronic event 
shape observables. We have not found a physically motivated modification of NNA that 
would alleviate this difficulty; therefore we investigated the ambiguity by comparing the 
higher-order coefficients obtained from restoring non-abelian contributions from both 
fermion and scalar loops.[] The comparison is shown in the last two columns of Tab. |I| 
in the V-scheme. To make a meaningful comparison in other schemes one would have to 
adjust the values of the strong coupling to take into account the difference between the 

^This comparison was suggested to us by Lance Dixon. 
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finite parts of the fermion (scalar) loop. (In the MS scheme C = —5/3 for quarks and C = 
—8/3 for scalars.) It is seen that the ratio scalars/quarks is very stable around 1.2, which 
is the asymptotic limit in large orders and corresponds to the ratio of the coefficients 
of in the expansion of the secondary quark distribution function. Consequently, the 
ambiguity in restoring the non-abelian contributions is at least 20%. However, in higher 
orders, the NNA prescription is probably less accurate than this anyway. The main 
point here is that radiative corrections to ctl have the same sign and are consistently 
large already in low orders. Keeping this in mind, we conform to restoring higher-order 
contributions from fermion loops in the following. 

The sum of {Po(^s)'^ contributions to all orders is conveniently written in terms of 
‘enhancement factors’ |^, measured relative to the leading-order contribution. They 
are dehned by 

OO 

(«,) = 1 + E , (94) 
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(95) 


In the following we dehne the sum to inhnity in the sense of a principal value Borel inte¬ 
gral, although truncating the series at its minimal term would be completely equivalent 
for practical purposes. For various values of as{Mz) we get, at Q = Mz, 

= 0.110 : MtPl = 1.59 mW = 1.92 ± 0.05. 


a. 


0.120 : 


= 1.68 


= 2.08 ±0.08 

0.130 : 


= 1.79 


= 2.23 ±0.12 


(96) 


The given uncertainties for the secondary quark contribution roughly coincide with the 
size of the minimal term in the series. The corresponding uncertainty for is small 
in comparison with the one for and is omitted. Interpolating between these values, 
one can get a rough idea of how changes with ag- 

Experience with similar calculations suggests that the approximation of resumming 
only (/dotts)"' contributions overestimates radiative corrections. A more realistic estimate 
would be expected to lie in between the resummed and the exact next-to-leading order 
(NLO) result. To see the magnitude of separate contributions, we take half of the Borel 
sum of (Poag)^ contributions starting at order n = 3 to estimate the size of higher-order 
corrections. With ag{Mz) = 0.118 we then get 

(7L/(Ttot = 0.0495 + (0.010 ± 0.010) ± 0.003 , (97) 

where the hrst number corresponds to the NLO result, the second bracket presents our 
estimate of further perturbative corrections, and the third gives the estimated uncer¬ 
tainty in summation of the perturbation theory. At Q = Mz this uncertainty is rather 
moderate in size, even though it corresponds to a 1/Q correction. This suggests that al¬ 
though non-perturbative corrections to ctl are much larger than to a tot, these corrections 
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Figure 9: Longitudinal fraction in the total e'''e“ cross section calculated with as{Mz) = 
0.118. Solid lines: leading order (LO), next-to-leading order (NLO) and resummation of 
all orders in corrected for the exact 0{a‘^) coefficient (NNA). Dashed line is the 

prediction of the Dokshitser-Webber model and dotted line is the average of NLO and 
NNA calculations. 

are still not large at Q = Mz- An exact 0(0;^) calculation would reduce the theoretical 
error considerably. Then, as determined from ul and atot would provide an interesting 
consistency check. 

Apart from measurements at a hxed energy, measuring the Q dependence of 
would be extremely interesting, especially at moderate energies. For other event shapes, 
it is known that a sizeable ‘hadronization correction’ must be added to second-order 
perturbation theory in order to reproduce the energy dependence. On the other hand, 
once higher-order corrections are computed, for example in an approximation such as 
‘naive non-abelianization’, such hadronization corrections must be reconsidered, since 
the sum of higher-order corrections can already produce a steeper energy dependence. 

To illustrate this point, we have plotted in Fig. || the energy dependence of the 
total longitudinal cross section, showing the leading-order (LO), next-to-leading order 
(NLO) and the resummed (NNA) results. We note the pronounced energy dependence 
of the resummed result towards lower energies. It is noteworthy that Monte Carlo 
parton showers do not produce such an energy dependence (see Fig. 5 of p|) and resum 
a different set of higher-order corrections. At the same time, since higher-orders in 
(/3oa<j)" are inseparable from 1/Q power corrections, it is plausible that a substantial 
part of the conventional hadronization correction effectively parametrizes these higher- 
order corrections. 

A similar conclusion can be drawn from a comparison of the resummed NNA re- 


35 




suit with the prediction obtained in the approach of Dokshitser and Webber (DW) [Q, 
compare the dotted and the dashed curves in Fig. The Dokshitser-Webber approach 
assumes universality of 1/Q power corrections and parametrizes them by an effective 
coupling. The dashed curve has been obtained following the prescription of with the 
effective coupling htted from the average 1 — T (T is thrust) and the relative coefficient 
oi 1/Q for 1 — T and determined in the massive gluon scheme. In both cases the 
perturbative expansions are trnncated at second order. It is seen that the energy de¬ 
pendence of the resnmmed pertnrbative results is not too different in NNA and the DW 
model. There is no conflict between the procedure of |]^ and the resnmmation presented 
here, if the phenomenological 1/Q correction effectively parametrizes the higher-order 
pertnrbative contribntions added in onr approach. If nniversality of power corrections 
holds, these perturbative corrections would also be universal, at least asymptotically in 
large orders. However, from the point of view presented here, the nniversality assnmp- 
tion is not reqnired, since higher-order corrections are in principle calculable for each 
observable. 


6 Conclusions 


In this paper we investigated power-suppressed corrections to fragmentation processes. 
Drawing on ideas derived from large-order perturbation theory and applied before with 
some success to deep-inelastic scattering |^, |2^, |2l| , we have modelled the x dependence 
of the leading 1/Q^ power corrections. The method we use is consistent with the QCD 
light-cone expansion of Ref. [§| and can be considered as a model of the x dependence 
of certain twist-fonr correlation fnnctions, based on the assnmption of ‘ultraviolet dom¬ 
inance’ of higher-twist matrix elements. The theoretical status of this assumption, as 
discussed in Sect. 2, is rather dnbious — eventnally experimental results should clarify 
its relevance, as in the case of deep-inelastic scattering. 

For the hrst time we were able to nnmerically analyse additional ambiguities that 
affect the method for observables where the hadronic hnal state is weighted or not fnlly 
inclusive. These additional ambiguities affect power corrections to event shapes in gen¬ 
eral [^] and also shed some light on the reliability of ‘massive gluon calculations’. In 
Sect. 3.5 we provided a parametrization of the x dependence of the leading (1/Q^) power 
corrections, obtained by extracting the generic dependences on x from the comparison 
between different implementations of the method. The parametrization depends on fonr 
constants which, although of order unity, should be determined experimentally. We 
should mention that we did not address kinematic effects dne to the hnite mass of the 
hnal-state hadron, which should be treated separately, as is done for target-mass effects 
in deep-inelastic scattering. 

From the theoretical point of view, fragmentation processes look like a promising 
place to pursne farther ontstanding qnestions related to the existence and interpretation 
of 1/Q power corrections in processes withont operator prodnct expansion, snch as Drell- 
Yan production or event shapes. We have seen that the light-cone expansion breaks down 
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when X ~ A/Q and quantities such as moments of the fragmentation cross section, which 
include the region of such small x, do not have a light-cone expansion. This situation is 
quite different from deep-inelastic scattering. In particular, the integrated longitudinal 
and transverse cross sections have a 1/Q correction. Fractional power corrections occur 
in fractional moments and in general the power behaviour of x-integrated quantities 
depends on the details of how the small-x region is weighted. 

Acknowledgements. We are grateful to Bryan Webber for discussions and for 
informing us on |]^ prior to publication. M. B. thanks Stan Brodsky and Lance Dixon for 
interesting discussions and the Institute for Nuclear Theory in Seattle for its hospitality 
while part of this work was done. V. B. is grateful to the Institute for Nuclear Theory 
in Seattle and the CERN Theory Group for their hospitality. 


References 


[ 1 ] 

[ 2 ] 

[3] 

[4] 

[5] 

[ 6 ] 

[7] 

[ 8 ] 
[9] 


DELPHI collaboration, P. Abreu et al, Phys. Lett. B 311 (1993) 408; contribution 
to ICHEP’96, Warsaw [paO 1-022]. 

OPAL collaboration, R. Akers et al, Z. Phys. C 68 (1995) 203. 

ALEPH collaboration, D. Buskulic et al, Phys. Lett. B 357 (1995) 487; ibid. B 364 
(1995) 247 (E). 

P. Nason and B.R. Webber, Nucl. Phys. B 421 (1994) 473. 


W. Furmanski and R. Petronzio, Phys. Lett. B 97 (1980) 437. 


P.J. Rijken and W.L. van Neerven, Phys. Lett. B 386 (1996) 422, |hep-ph/960443q ; 
Leiden University preprints INLO-PUB-09-96, |hep-ph/9609377 and INLO- 
PUB-10-96, |hep-ph/9609^ . 

G.T. Bodwin, S.J. Brodsky and G.P. Lepage, Phys. Rev. D 39 (1989) 3287. 

LI. Balitsky and V.M. Braun, Nucl. Phys. B 361 (1991) 93. 


H.D. Politzer, Nucl. Phys. B 172 (1980) 349; R.K. Ellis, W. Furmanski and 
R. Petronzio, Nucl. Phys. B 212 (1983) 29; R.L. Jaffe, Nucl. Phys. B 229 (1983) 
205. 


[10] A.V Manohar and M.B. Wise, Phys. Lett. B 344 (1995) 407, [hep-ph/9406392 . 

[11] B.R. Webber, Phys. Lett. B 339 (1994) 148, |hep-ph/940822^ . 

[12] Yu.L. Dokshitser and B.R. Webber, Phys. Lett. B 352 (1995) 451, |hep-ph/9504219. 


37 













[13] R. Akhoury and V.I. Zakharov, Phys. Lett. B 357 (1995) 646, |hep-ph/950424^ ; 
Nucl. Phys. B 465 (1996) 295, [hep-ph/9507253|; Phys. Rev. Lett. 76 (1996) 2238, 
[hep-ph/9512433| . 


[14] P. Nason and M.H. Seymonr, Nucl. Phys. B 454 (1995) 291, [hep-ph/9506317 . 

[15] G.P. Korchemsky and G. Sterman, contribntion to Moriond 1995: QGD and High 
Energy Hadronic Interactions [hep-ph/9505391 . 


[16] Short reviews are given in M. Beneke, |hep-ph/96092lR to appear in the Proceed¬ 
ings of IGHEP’96, Warsaw; V.M. Brann, [hep-ph/9610212| , to appear in the Pro¬ 
ceedings of DPF’96, Minneapolis; R. Akhonry and V.I. Zakharov, hep-ph/9610492| , 
to appear in the Proceedings of QGD’96, Montpellier. 

[17] M. Beneke, V.M. Brann and L. Magnea, SLAG preprint SLAC-PUB-7274, 
[hep-ph/9609266| , to appear in the Proceedings of QGD’96, Montpellier. 

[18] M. Dasgnpta and B.R. Webber, Gavendish preprint Cavendish-HEP-96/9, 
|hep-ph/9608394j, to appear in Nucl. Phys. B. 

[19] Yn.L. Dokshitser, G. Marchesini and B.R. Webber, Nucl. Phys. B 469 (1996) 93, 

[hep-ph/9512336| . 

[20] E. Stein et al, Phys. Lett. B 376 (1996) 177, |hep-ph/9601356| ; M. Manl et al, TU 
Mnnich preprint TUM/T39-96-29, |hep-ph/961230C| . 

[21] M. Dasgnpta and B.R. Webber, Phys. Lett. B 382 (1996) 273, [hep-ph/9604388| . 

[22] M. Beneke, V.M. Brann and V.I. Zakharov, Phys. Rev. Lett. 73 (1994) 3058, 
|hep-ph/9405304 . 

[23] M. Beneke and V.M. Brann, Phys. Lett. B 348 (1995) 513, [hep-ph/9411229| ; P. Ball, 


M. Beneke and V.M. Brann, Nucl. Phys. B 452 (1995) 563, [hep-ph/9502300 


[24] M. Nenbert, Phys. Rev. D 51 (1995) 5924, |hep-ph/9412265| . 


[25] G.P. Korchemsky and G. Sterman, Nucl. Phys. B 437 (1995) 415, |hep-ph/9411211 


[26] M. Beneke and V.M. Brann, Nucl. Phys. B 454 (1995) 253, [hep-ph/9506452 


[27] B.R. Webber, Lectnres at the Snmmer School on Hadronic Aspects of Gollider 
Physics, Znoz, 1994, |hep-ph/9411384| . 

S.J. Brodsky, G.P. Lepage and P.B. Mackenzie, Phys. Rev. D 28 (1983) 228. 


38 







































